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I. INTRODUCTION 

The symmetry associated with the so-called Laplace-Runge-Lenz (LRL) vector in New- 
tonian Kepler- Coulomb 2-body systems^, well-known for more than two centuriea^*^, is re- 
garded by many as 'accidental' or 'hidden'-"-. It has gained this adjective being associated 
originally only with 1/r potentials and because of the particular properties of bound states 
in these interactions - closed orbits and extra degeneracy of the energy states; and, also, 
because it had no apparent more profound source than its mere existence, in strong contrast 
with the common space-time Galilean or Lorentz-Poincare symmetries, which are clearly 
geometrical in nature. 

The realization, following Bacry et al^, Fradkin^° and Mukunda^\ that LRL-like vectors 
do exist and may be defined for all systems with rotational symmetry, changed this picture 
drastically. The generalized LRL vectors generate (together with the internal angular mo- 
mentum) o(4) or 0(3, 1) Lie-Poisson algebras just as in the Kepler-Coulomb case, extending 
the LRL symmetry into various systems, including, in particular, centrally symmetric ones 
with open orbits. 

These developments dissociate the two aspects which in the past seemed, following the 
Kepler- Coulomb case, to be inseparably characterizing the LRL symmetry - closed orbits 
and extra degeneracy on the one hand, and the existence of a constant LRL vector and 
the algebraic structure associated with it on the other hand. While the extra degeneracy 
continues to be peculiar to some interactions only, a LRL vector always accompanies internal 
rotational symmetry. 

Following the three seminal papers, generalized LRL symmetry was realized many times 



in the literature in various systems, including relativistic ones (see, e.g., Refs. llSHlSj), the 



MICZ system (electric charge + magnetic monopole)^"— , systems with spatial constant 
curvature^i^S, or other systems with velocity and/or time-dependent interactions besides 
central potentials^i (see also Ref. [s! for a more exhaustive list). Once the association between 
the LRL symmetry and a particular type of interaction has been dissolved, the symmetry 
may not be regarded any more as 'dynamical' and 'accidental'. Being an inseparable aspect 
of any rotationally symmetric system, the LRL symmetry must be of geometrical origin and 
nature. 

The first step towards identifying and understanding the physical origin of the LRL 



symmetry was done by Dahl, who discovered^^i^ some years ago for the classical (Kepler- 
Coulomb) LRL symmetry that its origin resides within the relativistic framework : Dahl 
has shown that the Newtonian LRL vector appears naturally in the computation of the 
Lorentz boost in the post-Newtonian approximation of electromagnetic or gravitational 2- 
body systems. Although it is only the Newtonian LRL vector that appears there this 
an essentially relativistic result, because it is of order 1/c^, vanishing in the full non- 
relativistic limit when the Lorentz boost becomes the Galilei boost. In essence, Dahl's 
result stems from the observation that the Newtonian centre-of-mass (CM) of an A^-body 
system X-^ = (^^ rriaXa) / (Yla "^a) ^^ ^^^ ^^ reference frame, which is not constant for 
relativistic systems, nevertheless its time varying part must be purely relativistic because 
Xn is constant in the non-relativistic limit. Explicitly computing this time- varying part of 
Xn, two independent solutions ensue, the difference of which is proportional to the LRL 
vector. 

The generalizability of the LRL symmetry a-la Bacry et al, Fradkin and Mukunda, on 
the one hand, and the realization of its relativistic origin by Dahl on the other hand, now 
call for unification. The purpose of the present paper is therefore to demonstrate, confirm 
and further establish the generality and universality of the Laplace- Runge-Lenz symmetry, 
bringing together the two approaches. 

We start (SecJII]) by introducing two propositions which encompass the essence of the 
LRL symmetry in general rotationally symmetric systems. These propositions provide 
a novel and very simple demonstration of the generality of the LRL symmetry together 
with a straight-forward tool to compute the Lie-Poisson brackets of the LRL vector in any 
rotationally-symmetric constellation. The generic properties and the consequences of the 
LRL symmetry are discussed in Sec JIIII together with the effect of the corresponding sym- 
metry transformation. 

The explicit construction of the LRL vector for general centrally-symmetric 2-body sys- 
tems starts to be discusses in Sec JIVI The method for constructing the LRL vectors for 
centrally-symmetric systems with open orbits suggested in the pastp^s^i deemed these vectors 
to be only piece-wise conserved, changing directions non-continuously in turning points^i^i-JI^ 
thus reducing substantially any interest in their applicability and usefulness. In fact, almost 
all the known applications of generalized LRL vectors are in systems for which the vector is 
constant, corresponding to closed orbits. Since the LRL vectors that are obtained via Dahl's 



(relativistic) procedure (see Refj22l and Sec JVIl below) are constant, while even for simple 
relativistic systems the orbits are opert^^i^^, this situation is not satisfactory. To remedy the 
situation, we provide a definition for constant LRL vectors in general centrally symmetric 
systems that is valid also for open orbits, together with the explicit computation of its self 
PB. This method is then applied in Sec|V]to two relatively simple relativistic systems with 
open orbits, the relativistic Coulomb system and post-Newtonian electromagnetic or gravi- 
tational system. In both cases the LRL vector is explicitly constructed, via two alternative 
ways. 

Finally, we refer in Sec JVIl to the relativistic origin of the LRL symmetry in general 
centrally symmetric systems, showing that the Newtonian LRL vector found in the preceding 
sections is indeed derivable from the Lorentz boost of the corresponding post-Newtonian 
extensions. The article then concludes with a discussion and summation of the generic 
properties of the LRL symmetry. 

II. GENERALIZATION OF THE LAPLACE-RUNGE-LENZ SYMMETRY 

To elucidate now the way the LRL symmetry appears, in a natural way, in all rotationally 
symmetric systems, and to produce a tool which will greatly simplify the discussion of the 
symmetry in general systems, let us introduce two propositions which sum up the main 
properties of the LRL symmetry, generalizing the properties of the symmetry in the classical 
Kepler-Coulomb case. 

Rotational symmetry is assumed to be characterized by the existence of constant internal 
angular-momentum vector C. (the term internal refers here and in the following to any 
dynamical quantity which depends only on the relative coordinates of the particles and 
their relative motion and is invariant under uniform global translations; in the case of a 
single particle with a fixed centre of force, all internal quantities are defined relative to the 
centre of force.). It should be noted that for non-centrally symmetric interactions, say with 
spin or velocity dependence, d, is not just the orbital angular momentum with terms like 
r X p but includes also extra terms (as in Eq. flT^ below). 

The systems under consideration are also assumed to be endowed with Lie-Poisson brack- 
ets (PB) {., .}. For the general discussion in the present Section and Sec. [TTllit suffices that 
these PB satisfy the general requirements from Lie-Poisson brackets^^, and no canonical- 



simplectic structure needs to be assumed (such a structure, with the standard definition for 
the PB, will only apply to particular examples later on). 

The identification of C. as the generator of internal spatial rotations, relative to the centre 
of mass or centre of force, is then incorporated in the requirement for the existence of 
rotational PB 

{K\ f} = {t, K^] = e'^^K^ , (1) 

for any internal vector K. The anti-symmetry of the self PB {K^,K^} implies, for any 
internal vector K, the existence of another internal vector A such that 

{K\K^} =6'^''A'' (2) 

— * 

From Eq.(|2]) it follows, using the vector property ([1]) for K, that 

lK\i-K\ = {K\f}K^ + {K\K^}f = e'^'^A^f (3) 

The product i- K is therefore ^-invariant, in the sense that < K\ i ■ K > = 0, iS the vectors 
A and i are parallel, say as A = ai. We may therefore conclude that 

proposition 1 The self PB of an internal vector K are of the form 

{K\K^]=ae'^^t (4) 

with a some scalar, iff the product £ ■ K is K -invariant, <K^,i-K> =0. 

In 2-body systems any scalar observable constant of the motion must be functionally 
dependent on H and i^. Hence, if K^ is a constant of the motion then it must be a 
function of H and v^, say K'^ = F(H,i'^). For more general systems, constant scalar 
observables need not be functionally dependent on H and £^ only. However, we may still 
consider those vectors K for which the product i ■ K is i^-invariant, {K,i ■ K} = 0, and 
their magnitude K^ depends, besides H and £^, also upon other observables which are i— 
and i^— invariant. These observables will henceforth be generically denoted as A (thus 
li,A\ = 0, |i?,^| = 0), so that we may write now i?^ = F{H,f,A). Then we have 
proposition 2 Let K he an internal vector such that : 

1. The product i ■ K is K -invariant 



2. K^ = F{H,e,A) 
Then the self PB of K satisfies 



d K' 



{K\K^] = -^^e^H' (5) 

Proof The PB < i^*, K^ \ may be computed either as 

iK\K'^\=2{K\K^]K^ = 2ae'^^KH'' (6) 



{k\K^} = {K\F{HA\A)] = —^ {K\f} = ^^r^f'PK' (7) 



or as 

Comparing the last two equations and taking into account the non-parallehsm of i and K, 
Eq.® follows. QED 

To illustrate the applicability of these propositions, here are few examples : 

1. First, for Newtonian 2-body Kepler- Coulomb systems with the Hamiltonian 

«-f + - (8) 

2/x r 

{f = Xi — X2 is the relative coordinate, p the corresponding momentum and /x the 

Newtonian reduced mass) and internal angular momentum i = fx p, the LRL vector 

is commonly defined as 

K = pxl+—r (9) 

r 

Substituting the magnitude K^ = 2fiHi'^ + fi^K^ in Eq.(|5]) then yields the well-known 
result 



{K\K^} = -2iJiHe'^^t (10) 

2. In systems with a magnetic monopole and general central interaction the equation of 

motion is 

m—- = —V X r U [r) r (11) 

at r-^ r 

with conserved energy and angular momentum 

X —* (y. 

E = -mv^ + U (r) , i = mfxv r (12) 

2 r 



In particular, for a MICZ^"^ system with modified Coulomb interaction 

CM-F + £^ (13) 

the LRL vector takes the particularly simple form 

^ -» rriK, ^ , . 

K = mv X i -\ r (14) 

r 

Here K ■ i = —maK : non-zero, but being composed only of constants of the system 

it is i^-invariant. Substituting the magnitude K^ = 2mE{£'^ — a^) + m'^K,'^ in Eq.(l5]) 

then yields immediately, without any tedious computations, 

{K\K^} = -2mEe''''i'' (15) 

In conclusion, we note that : 

1. As is evident from its derivation, Eq.(|5]) is independent of any particular form of the 
interaction, except for the requirement for rotational symmetry. 

2. Eq.([5]) provides a simple and straight-forward tool to compute the self PB of K. Its 
usefulness will become evident in the various applications in the following. In partic- 
ular, comparing with Fradkin's^ computations for LRL vectors in general centrally 
symmetric systems, Eq.(l5]) will prove to save need to launch on very tedious calcula- 
tions. 

3. Moreover, Eq.Q doesn't require any particular recipe for the computation of the PB. 
It is therefore also suitable for systems which lack canonical or phase-space structure, 
such as in some relativistic act ion- at- a- distance systems. 

Any internal vector observable K that satisfies these conditions may be regarded as a 
generalized Laplace- Rung e-Lenz vector. We now proceed for further applications of these 
propositions, to establish and confirm the universality of the LRL symmetry. 

III. ALGEBRAIC PROPERTIES AND CONSEQUENCES OF THE 
GENERALIZED LAPLACE-RUNGE-LENZ SYMMETRY 

From the discussion of Sec. [TTl it follows that any constant vector which is not parallel to 
£ could serve, at least in principle, as a LRL vector. This becomes evident in the following 



way: Let Uo be a constant unit vector perpendicular to L By Eq.([5]), any vector whose 
square is independent of I has vanishing self PB. Therefore, any such constant unit vector 
is a generalized LRL vector with vanishing self PB, 

Mo ■ f = , Ua-Uo=l => {ui, Ui} = , (16) 

and may be used to generate arbitrary generalized LRL vectors via the relation 

K = f{E, f, A)uo + g{E, f, A)ix u^ + ^^^^^/ (17) 

The coefficients f(E,i'^,A), g{E,i'^,A) and h(E,A) may be arbitrary functions of their 
arguments, where A stands for possible constant scalar observables which are i- and K- 
invariant. K ■ i = h(E,A) cannot depend on i^ since it must be ^-invariant. 

The LRL vectors constructed in Eq.l lTTll are very general, and in principle, with an 
appropriate choice of the coefficients one may get almost any desired result. However, the 
physical meaning of such a construction would be obscure. We therefore seek now to limit 
this generality by introducing some physical content. 

By properly choosing the coefficients in Eq. flT7|) and taking Eq.([5]) into account, a LRL 
vector A with self PB 

{A\A^}=r]e'^''i'' = r]t^ (18) 

may always be created with the coefficient t] being either +1 or —1 (the critical limiting 
value r] = needs not be considered separately). Adding the PB of rotation, 

[A\ f] = {f , A^} = e'^^A^ , (19) 

the vectors A and £ generate together classical Lie-Poisson algebras o(4) or o(3, 1) according 
to wether rj = +1 or —1, respectively. Any vector A satisfying Eq.( TT8|) may be regarded as 
a canonical Laplace- Rung e-Lenz vector. 

In quantum systems, the Casimir operators determine the quantum state of the system. 
In classical systems, the Casimir invariants of the corresponding Lie-Poisson algebra provide 
information regarding the physical state of the system. With the PB ( TTSll and (fT9|) . the two 
Casimir invariants are 

Ci = Ci {H, A) = r]A^ + ? (20a) 

C2 = C2{H,A) = A-e (20b) 

8 



For many 2-body systems, in particular centrally-symmetric ones, the motion is in a plane 
perpendicular to L Then the LRL vectors, and in particular the canonical one A, may be 
chosen in the plane of motion so that C2 = A ■ £ = 0. Other systems, such as the MlCZr^ii^ 
systems discussed in example B in Sec. [Ill [Eqs. flllH14p ] or generalizations thereof such as 
those discussed by Iwai and Katayama^, may be characterized by non-zero €2- However, 
even in the latter case it follows from Eq. (fT71) that the component of K parallel to £ may be 
chosen at will. Therefore, if C2 7^ 0, we may always use the component of A perpendicular 
to £ instead of A. In other words, it is always possible to assume C2 = without loss of 
generality, and this assumption will be held in the following. 

The sign coefficient r] is not arbitrary, but reflects the energetical state of the system as 
for classical Kepler- Coulomb systems. For a given value of the total energy E, a system may 
be either bound or un-bound. The o(4) case (r/ = +1) corresponds to bound states : From 
Eq.f l20ap it follows that Ci is necessarily non- negative, with (? bounded from above, 

e < Ci (21) 

The existence of an upper limit for i (for a given value of E) is characteristic of bound 
states. In particular, in the case of circular motion rotational invariance requires A = 0. 
Then, with given total energy E, i achieves its maximal valued. Denoting this maximal 
value as irnaxiE) it determines Ci as^ 

C,{E) = el,^{E) (22) 

The opposite case (r/ = —1), with the o(3, 1) algebra, corresponds to un-bound states: Since 
y4^ = £^ — Ci > there is no limitation on i, but Ci must be negative in order to allow 
situations like head-on collisions in which i = 0. Still, except for the sign, the functional 
dependence of Ci on E is the same as in Eq. fl22|) . 

Another aspect of the general construction of the LRL vector in Eq. lfTT]) is that its 
direction may also be chosen at will. It is evident that the PB of A in Eqs. f lTSj) and 
(II 9p do not change if A is rotated perpendicular to i. Therefore, at least as far as the 
algebraic relations are concerned, the LRL symmetry does not distinguish any particular 
direction for the LRL vector. 

There is, however, a way to determine a preferred direction for the LRL vector : From the 
post-Newtonian derivation by DahP- emerges the classical LRL vector ([9]), directed towards 



the perihelion. Later, in Sec JVIt it is verified that the LRL vector that emerges from a 
similar derivation for general centrally-symmetric systems also points in the direction of 
closest approach ('generalized perihelia'). 

We conclude the present section with a discussion of the effect of the LRL symmetry on its 
own algebra. Symmetries in classical dynamics were firstly manifested in terms of mappings 
of the configuration- or phase-space. This has certainly been the case with the common 
space-time symmetries (translations, rotations, dilatations, etc.). In contrast, with the LRL 
symmetry the focus has always been on the constant LRL vector and the corresponding 
algebraic structure that it generates together with the internal angular momentum. The 
transformations which constitute the LRL symmetry group for classical Kepler- Coulomb 
systems were discussed^"— only many years after the discovery of the algebra that dominates 
the symmetry. 

In the foregoing discussion we considered the generic properties of the LRL symmetry, 
so generic even to the extent of not assuming any particular structure for the configuration- 
or phase-space. But even in the absence of such structure, the algebra can provide an 
insight into the nature of these transformations. While £ generates internal rotations, the 
internal transformations that correspond to the LRL symmetry are generated by A. Let n 
be some constant unit vector and x be a real dimensionless parameter. The infinitesimal 
transformations generated by A relative to the direction of n via the generator 5Q = n- A5x 
satisfy the equations 



51 



5g\ =nyi A5x 



5A=\A 



,5Q\ = ijn X i6x 



(23) 



The orbits of these transformations within the algebra, in terms of A and i, are given by 
Rodrigues-like formulae 



r/ = +1 : 



?7 = — 1 



£{x) = cos X ■ £{x = 0) + (1 ~ cos x) £{x = 0) ■ n n + sin xn x A{x = 0) 
A{x) = cos X ■ A{x = 0) + (1 — cos x) A{x = 0) ■ n n + sin xn x i{x = (X24a) 

£{x) = cosh X ■ £{x = 0) — (cosh x ~ 1) ^(x = 0) ■ n n + sinh xn x A{x = 0) 
A{x) = cosh X ■ A{x = 0) — (cosh x — 1) A{x = 0) ■ n n — sinh xn x i[x =(2)|b) 



10 



These transformations, relating orbits with equal energy but different angular momentum, 
may be regarded as deforming transformations. They are the non-quantum analog of the 
shift operators generated by the LRL vector for the hydrogen atom wave function^"—, 
extending to arbitrary centrally symmetric potentials. 

A-generated transformations may change the direction of L For centrally-symmetric 2- 
body systems this implies changing the orientation of of the plane of motion, which is merely 
a geometrical transformation rather than physical. In such systems we may always choose 
the vectors A and n perpendicular to d, thus restricting the algebra, without loss of generality, 
to the subalgebra defined by C2 = A ■ i = 0. The vector i may then change its magnitude 
but not the direction, keeping the plane of motion intact. The physical system is then fully 
characterized by the first Casimir invariant Ci, while different states or configurations are 
characterized by i and A. 

In the case of bound centrally symmetric systems it is convenient to start, at x = 0, 
at the state of circular motion for which A = 0. We may be interested in particular to 
maintain the direction of the angular momentum unchanged, which implies that n must be 
perpendicular to i. Then, using i{x = 0) = iuia.x{E)£, we obtain from Eq. fl24al) 

A{X)= ^max(^) sin xnxi = tan xn x i= ^ il^^^{E) - P n x I (25) 



IV. EXPLICIT CONSTRUCTION AND COMPUTATION OF 
GENERALIZED LAPLACE-RUNGE-LENZ VECTORS 

Fradkin^^, and later also Perea^ and Yoshidai^, suggested and discussed a general method 
for explicitly constructing the generalized LRL vectors for general centrally-symmetric sys- 
tems with open orbits. This method of construction deemed these vectors to be only 'piece- 
wise' conserved, changing directions non-continuously in turning points^i^"— , thus reducing 
substantially any interest in their applicability and usefulness. Consequently, almost all 
the known applications of generalized LRL vectors are in systems for which the vector is 
constant, corresponding to closed orbits. However, taking into account the fact that the 



LRL vectors that are obtained via Dahl's (relativistic) procedure (see Refl22l and Sec JVII be 



low) are constant, together with the fact that even for simple relativistic systems the orbits 

11 



are open^^i^^i^, and the results of the Sees. HTl and IIIII implying that constant LRL vectors 
should exist and be definable for any rotationally symmetric system, this situation calls for 
a remedy. 

The purpose of the present Section is therefore to provide an improved definition of the 
generalized LRL vector, insuring its constancy for any relevant system. Applying Eq.dS]), an 
expression for the self PB of the LRL vector is readily obtained. 

In the following we consider 2-body centrally symmetric systems, with canonical internal 
variables (r, p) and conserved angular momentum i = fxp. The construction of LRL vectors 
in general systems commonly starts^'^^'^i with choosing an arbitrary constant unit vector 
Mo, in any direction perpendicular to i, as in Eg. (11 71) . As we have shown, any such vector 
satisfies Eq. flT6l) and may be regarded as a LRL vector. With the motion being confined to 
the plane perpendicular to i, let 6 be the azimuthal angle defined counter-clockwise from Uo 
to f. In terms of Uo and 6, the unit vector f = f/r may be represented as 

f = cos 6uo + sin 6i x Uo = U{6)uo (26) 

where d, = i/i and 

u{e) = cose + sin eix (27) 

is the rotation operator which rotates vectors in the plane of motion counter-clockwise with 
angle 6. Inverting Eq. fl2B]) . it is possible to express Uo in terms of f and 6, 

Uo{f, e) = U{-e)f = COS 9r - sin Oi x f (28) 

Taking the constant unit vector Uo as the direction of the desired generalized LRL vector, 
and using arbitrarily any scalar function K{E, i'^) for its magnitude, we obtain the vector 

K = K{E,f) ■ Uo = KU{-9)f = K cos 9f - KsinOi x f (29) 

which satisfies all the requirements from a generalized LRL vector. Using the identity 

^ i'^f— r'^px £ 
i X r = — — 



r ■ p 

and the notation pr = "r ■ p, it may be brought to a more familiar form 

^ Ksine ^ / Kisine\ , , 

K=—„ pxi\\Kcosd f (30) 

ipr V "^Vr J 

12 



which resembles the classical LRL vector ([9]). 

Though similar to the construction of Fradkin^ and Peres^i (shown to be equivalent by 
Yoshida^^), this construction differs in one essential aspect - the scalar coefficients oi px i 
and f in Eq.f l30p are functions of the polar canonical variables (r, 6]Pr,pe = ^), while in the 
construction of Fradkin and Peres the corresponding coefficients are only functions of r, 
after having integrated the equations of the orbit for the particular orbit. Our construction 
need not take into account the equations of motion, and it doesn't suffer, therefore, from 
the drawbacks of the Fradkin-Peres construction - namely, being only piece-wise conserved. 
In fact, the construction in Eq. fl30|) is just a complicated way to write the constant vector 
K = K{E,i'^) ■ Uo ( 129|) which, by virtue of Eq.([5]), is already recognized as a LRL vector. 

So far, the direction of the unit vector Uo in the plane of motion is completely arbitrary. 
In order to fix it, we now apply the equations of motion. 

For an arbitrary central potential the orbit may always be represented as r = r (6*, £", C). 
This is a periodic function in ^, with period G being some general function depending on the 
parameters of the orbit, 6 (-E, €). Expressing r as a function of 9 via r (6*, E, i), the coefficient 
oi px i in Eq. fl30|) is then in general a varying function of 6, constant only for 1/r potentials, 
as in Eq.(|9]). Although the vector (130|) is certainly regular, by its construction, the separate 
coefficients there may be singular because Pr vanishes at the turning points of the orbit. 
In order to make the coefficients regular we require that sin 6' also vanish there, which is 
achieved by choosing Uo so that it points towards a turning point. The speciality of the 1/r 
potentials is that only in this case = 27r so sin 6' may be made to vanish simultaneously 
at all the turning points; for all other potentials sin 6 may be made to vanish at one turning 
point, but there will be other turning points (if the system is bound) for which {sm6)/pr is 
singular. 

To obtain explicit expressions in Eq.( l30|) for particular systems, we start from the fact 
that in any non-circular configuration of the systems in question there is at least one point of 
minimal approach. For unbound systems there is just one point like this, which is the only 
turning point. For bound systems, except for 1/r potentials, there are multiple points of 
minimal approach, so we choose arbitrarily one of them. Let the chosen minimal approach 
point be at r = fv, directed from the centre-of-mass. 

Due to the central symmetry it is convenient to use polar coordinates, with an Hamilto- 
nian of the most general form, H = H {r,Pr,Pe) and pe = i is a. constant since dH/dO = 0. 

13 



Assuming symmetry under spatial reflections and time reversal, H must be an even func- 
tion in pr and pg. Let us choose Uo = rm-, the unit vector along fvn, so that 6* = there, pr 
vanishes at all the turning points since r = dH/dpr is odd in p,., and r^ is determined as a 
function of E and £^ as the smallest positive root of the equation 

H{r^,Pr = 0,i) = E ^ rm = rm{E,e) (31) 

A-priori, the magnitude of K may be chosen at will. To fix it with a vector K that resembles 
the classical LRL vector (jH]) the most we introduce the condition that the coefficient oi px £ 
in Eq. (15U]) be equal to 1 at r = fVi, 

■Ksine\ ^ ^32) 



SO that K is determined by 



r=rm,e=0 

Then, using Hamilton's equations from which follows that 

dpr dH/dr 

'de^^dH/di' 

the resultant LRL vector is 



^■(^•^^'^<^).......„ (^^' 



(34) 



These results apply to any centrally-symmetric system. For more specific results, let us 
consider Newtonian systems with Hamiltonian 

H = ^ + U{r) (36) 

The equation for the minimal distance is 

(p?)_„ = 2/.[i?-f/(r„)]--^ = (37) 

' m 

SO that the magnitude of the generalized LRL vector becomes 

K{E,e)= — -^JirJU'{rr^) (38) 

14 



and the resultant LRL vector is then 

02 



K-- 



Tr 



- firm U' (r. 



HTj^U \Trn) 



COS 6r — sin 6i x r 



(39) 



For 1/r potentials the classical LRL vector (Q is evidently obtained in Eg. (15^ . We 
also note that for a general potential U{r), in the case of circular motion for which pr = 
everywhere, Vm [E, £'^^^{E)] = To is the radius of the orbit. Then also K vanishes identically, 
K [E,f^.^^{E)] = 0. Otherwise, since r > r„, fsT > 0. 

The self PB of the vector K in Eg. (135!) may be computed by application of Eg. ([5]) and 
using Egs. flHTl) and fl5^ . In particular, for centrally symmetric Newtonian systems we get 
from Egs. ^ and ([MD 

4tti ' 



OK 

dW). 



- /i [rm U' (rm)] ^fiE - ji [rm U (r^)] 



2rrr,^K 



2K 



so that the self PB of the vector K (Eg. (!39|) ) are 



^ ' ^" din 



2^E-f,[rlUirJ]"]e^^H' 



(40) 



It is interesting to note that Eg.l lTOl) is obtained not only for 1/r potentials, but for all 
1/r + 1/r^ potentials. This is due to the fact that the 1/r^ term causes the conic section, 
which is the orbit for the 1/r term alone, to rotate in constant angular rate but leaves the 
shape of the conic section intaci^i'^. 

Another way to construct the generalized LRL vectors, which is useful for the post- 
Newtonian computation in Sec lVIl starts with the interaction-free part px i of the classical 
Kepler- Coulomb LRL vector (Q. With the Hamiltonian ( l36l) its time derivative is 



d 
~dt 



Let W{r,e,pr 



p X 



,_jji ^j,^f?y^£— _jj' (^^^ Ir'^p— (r ■ p) r] 



^rU' [r] 



df 



d 



PrU' {r)f= — [firlf' (r) r\ — Pr [rU (r)]" f 



dt dt 

be a vector observable which satisfies 

dW 



dt 



Pr [rU (r)] r 



with the initial condition W {9 = Q) = Q. Then the vector 

K = px £ — urU' {r)f+ W = [p'^ — ixrU' (r)] f — [p- f)p + W 



(41) 



(42) 



(43) 
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is clearly a constant of the motion. For 1/r potentials the RHS of Eq. (H2|) vanishes identically, 
and the vector in Eq. fH3|) reduces to the classical LRL vector (|9]) with W = 0. Otherwise, 
for a general potential U (r), the two vectors in Eqs. (15^ and (H5]l coincide at r = fvi; 
and, being both constant, are necessarily identical. This establishes the existence of W, and 
comparing the vectors in Eqs. (H5]) and (15UI) . with K in the latter given by Eq.( 15H]) . may be 
used to obtain an explicit expression for W. 

It is instructive to verify directly that W is indeed regular in the neighbourhood of fm- 
Converting the time derivative in Eq. fj42|) into an angular derivative using the equation of 
motion for 6, 






lir^ 



we obtain 



dW 



fir'^Pr 



rU (r)] r 



Pr vanishes at the turning points, but its angular derivative 

dpr i fir'^U' (r) 

'de^ r 1 ' 



(44) 



(45) 



does not (for non-circular configurations). Then pr = 0{9) near f = fm and consequently 
W{e) = O (6^) there. 

We conclude this section with a recipe for the PB of the generalized LRL vector K = Kuo 
with any desired observable J-". Towards this end, it is most convenient to use Eq. fl28|) for 
the representation of the unit vector Uo as an observable. The computation of derivatives of 
functions of 9 is performed using the geometrical relation 



de 



X r] ■ df 



(ixfY 



(46) 



and we obtain, after some algebra, for the PB of any observable J-" with Uo, 

sin 6 



{T, Uo} 



P 



or or 



(47) 



Hence, the PB of J-' with K given by Eq. (129|) are 

K KsinO 



{^,k} 



""U^.m + ^l^.n 



dH 



d{P 



K 



£3 



dF 



dJ^\ 



r^—+Pr^7^ 

or op J 



(48) 
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V. THE GENERALIZED LAPLACE-RUNGE-LENZ SYMMETRY IN 
SOME EXEMPLARY SYSTEMS 

The LRL symmetry for general Newtonian centrally symmetric systems was discussed in 
Sec. IIVI The relativistic origin of the LRL symmetry, first realized by Dahl^^ for Coulomb- 
Kepler systems and generalized to arbitrary centrally-symmetric systems in the following 
(Sec JVI|) . leads us to consider its appearance in relativistic systems with particular interest. 
Relativistic systems, even simple ones like the relativistic Coulomb system or the post- 
Newtonian systems, are characteristically endowed with bound states with open orbits. Still, 
the LRL vector that emerges from Dahl's procedure is constant. Thus it was important to 
insure in Sec. HV] that we are equipped with a valid definition for constant LRL vectors. 

So far, when the LRL vector was considered in a system with open orbits, the attitude 
was to start with the classical LRL vector (Q and to follow its rotation—"—. The constant 
LRL vectors computed by the method developed above certainly reduce to the classical LRL 
vector in the limit of closed orbits, and otherwise, already contain, built in, the data about 
the rotation of the orbit. We now derive and discuss, as a demonstration and application of 
the results and methods developed above, the generahzed LRL vector in these two relativistic 
systems. 

A. Laplace-Runge-Lenz symmetry in relativistic Coulomb systems 

Let us consider a relativistic Coulomb system with one of the particles having infinite 
mass and located at rest at the centre-of-mass. The dynamics of the other particle with 
mass m are determined by the Hamiltonian (in the present Subsection the convention c = 1 
is used) 



H = v/p2 + m2 + - = Jp2 + ^ + m^ + - (49) 

This is a well-known text-book problen>2^ characterized by irregular orbits, including open 
orbits for bound states^^. So far, only rotating LRL vectors were considered^"^, attached 
to the axis of the orbit. In the following we construct the constant LRL vector. 

In a configuration with given energy H = E and internal angular momentum pg = i, the 
squared radial momentum is isolated as 



Pt = {E - ^Y -m' -^-^ = E' -m' -^-^ - ^^ , (50) 
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The type of the orbit depends on the relative values of E and M, i and k. The major 
relativistic effect is the reduction of the centrifugal barrier from the Newtonian value i'^/r'^ 
to (£^ — K^) /r^. Hi > \k,\, so the centrifugal barrier is still maintained, the orbits are 
precessing conic sections. This is a well-known text-book result^. If the centrifugal barrier 
disappears {i = \k,\) or even reverses {i < \k,\, becoming kind of "centrifugal propeller") the 
orbits become irregular—. 

The LRL vector is determined by the turning points, which are determined by the equa- 
tion Pr = 0. In those orbits in which there is a distance of closest approach (otherwise simply 

r > 0), it is given by 

_Ek+ ^{E^-m^)P + m^K^ 
E^ — m^ 

Once the direction of closest approach is fixed (in some arbitrary direction), the LRL vector 
K is determined by Eq. (l35l) with the magnitude 

K {E, f) = -i (Hyl^) = V(i?2-m2)£2 + ^2^2 ^ (52) 

computed from the Hamiltonian fH9|) and using Eq.( l5T|) . Its self PB are then found, applying 
Eq.05]) and using Eq. (l52|) . to be 

{K\ K^] = - {E^ - m^) e'^^t (53) 

The type of the algebra is therefore determined solely by the energetic state of the system; 
although the forms of the orbits depend crucially also on the value of i'^ — n"^, the algebra 
is independent of the centrifugal condition. The LRL transformations i^^ may therefore 
take the system across the critical point oi i = \k,\ without difficulty. The magnitude of the 
canonical LRL vector satisfies 

^'iE^^')=T^?-^-V^'^ (54) 



1^2 



m 



with coefficient 1] = — sign (ii^^ _.^2^^ ^^ complete agreement with Eq.f jTSj) . and the first 
Casimir invariant being 



m^K^ 



^i(^) = TW-2 ^ (55) 

\E^ — m^l 

It is easily verified that for bound states Eq. (!22|) is satisfied, since for circular motion 

Ca.(i?) = -J== (56) 
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and otherwise d < iraa.x{E). Since i = r xp and the motion is always in a plane perpendicular 
to i, K ■ i = and the second Casimir invariant C2{E) vanishes. 

The generic form of the LRL vector is given by Eg. (1301) . Since the orbits are open, some 
dependence of K on the azimuthal angle 6 is unavoidable even when details of the orbit are 
taken into account. More concrete expressions for the LRL vector depend on the type of 
orbit. Using the relation 

which follows from Hamilton's equations, the equation for the orbit is obtained from Eq.f lSUp . 

I drV / fi;2\ 1 2kE _E'^-m^ 

Assuming that i'^ > k'^, so that a centrifugal barrier does exist, the solution for the orbit is 



1 Ek K{E,P) I kW 

with K {Eji"^), the magnitude of the LRL vector, given by Eq. (!52|) . Defining the angle 



k2 



^^dl--9 



the orbit is a conic section which is fixed in the r — (p plane, but rotating in the r — 9 
plane. To account for the rotation, lei ip ^ 9 — ip. Using the rotation operator identity 
U{a + 13)= U{a)U{/3) and Eq. jl^U]) . the LRL vector may be written as 

K = U{-9)rrn = U{-^ - ip)r^ = U{-^)U{-^)f^ 

Then, using the relations 

f - K^ 

Kcos(p = \- kE (60a) 



k"^ (?■ dr 



Ksin^ = yi - ^-,^^ = VW^^'Pr (60b) 

which follow from Eq. (!59|) together with Eq. (!57|) . the vector 

K' = U{-if)r^ = — - — -p X e+ [ K cos ip r = 

ipr V '^Pr ) 



K,^ ^ ^ I —K^—tyJt^ — K 



e-K^- (iy/W 



2 _ / /72 _ ^2 _ /) //72 _ ^2 



l--^px£+( -^ + KE]f (61) 
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is identified as the constant LRL vector in the fictitious r — ip plane, but it rotates together 
with the conic section in the physical r — 9 plane. K' coincides with the result of Yoshida^, 
obtained in a much more complicated way. Thus, finally, the LRL vector may be written, 
for (? > K^, as 



K = U{-i)) 



K 



2 



e -t^- l^P - H? 



1 - ^P X £ + I + kE ) f 



(62) 



For the critical value £ = |k| for which the centrifugal barrier disappears, Eq. (!58|) becomes 

( I drV 2E _ E^- m^ 

with the solution 

2kE 

^ ~ E^-m^- EW ^ ^ 

The LRL vector ( l62l) then reduces to 

K = U{-e)KEf (65) 

In the anti-centrifugal case £ < |fi:| the solution for the orbit becomes 

(66) 



Ek-Kcos\i\ A/€-ie 



Expressing 6 in terms of r and substituting in the generic LRL vector (I5U]) . will then yield 
the LRL vector in terms of r. 

B. The generalized Laplace-Runge-Lenz symmetry in post-Newtonian 
Kepler-Coulomb systems 

As a second example we derive in the following the generalized LRL vector for a 2-particle 
system in the post-Newtonian approximation either for electromagnetic (Darwin^^i^) or 
gravitational (Einstein-Infeld-Hoffman^^i^) interactions. As for the former case, LRL vectors 
were considered in these systems only as rotating classical LRL vectors (see, e.g.,—). The 
only time that constant LRL vectors were considered for these systems was, to the author's 
best knowledge, by Argiieso and Sanz^^, who computed the LRL vector by quite cumbersome 
means. Our derivation, applying the methods developed in the preceding sections, is, to our 
belief, more elegant and illuminating. It is also used to show a different way of application. 
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Consider a 2-particle system with masses mi,m2 and possible electrical charges 61,62- 
The post-Newtonian Hamiltonian for the two interactions in the CM system with internal 
polar canonical variables is given by2^ 



H 



p^ 



p 



K 

+ - + 



{2 + a)pl + {l + a) 



2yU 8z/'^c^ r 2mim2C^r 
with p"^ = pI + pl/r"^, Mo = rrii + 1712, fi = mim2/Mo, and 

1 _ 1 1 _ 1 8 



+ 



an, 



6Moc\ 






(67) 



r/i'f 1712 /^^ MolJ? 
K = 6162 and a = or K = —Gmim2 and a = 3Mo//i for electromagnetic or gravitational 

systems, respectively. 

Instead of proceeding as before, namely computing r^ from Eq.( l3T]) and then computing 

K via Eq. fl35|) . we use a slightly alternative way. With constant H = E' and pg = i, pi is 

explicitly deduced from the hamiltonian (l67Il 



^?-^--^K-S)? 



1^1 — a) /i 
Mo 



2kE' 



(3 + Sa) /i 



3Mo 



+ 



(68) 



From Hamilton's equations 

dr OH 

dt dpr 

de _ dH 

dt dpe 



p^ {2 + a)K 



+ 



/i 2v^c^ mim2C^r 



'l + a)K 



Pr 



p^ 



+ 



/i 2z/'^c^ mim2C^r 



it is easy to see that Pr may be written as 



du 



where, correct to terms up to order 1/c^, 



K 



_ 1 



(69) 



(70) 



Substituting ( 169|1 and ( ITOll in Eq.( l68|) and keeping terms up to order 1/e^ then yields the 
equation 



;i-*)-M'|) +(«-".)^ = b^ 



(71) 



where 



Mo 



B^ 



1+1 



1 + 



^3iV^ 
(So; — 3) /i 



/iK 



'1 — a) /i 



Mn 



kE' 



3Mo 






+ 1 



3fi\ E 



/2 



Mo / £2c2 



+ 



1+ 1 + 



5afi\ 2k, 



3MoJ Pc^ 



flK, 



1 + 



5a; /i\ K 



3MoJ 2£2c2 
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The solution, with the condition that ^ = at a perihehon, is 



u = Uo + Bcos[{l -5)9] 



(72) 



which immediately yields 



Pr = -i 



„du 



'l-6)iBsm[{l-6)9] 



(73) 



and the magnitude K {E', f) [Eq.dSS]) 



K 



dpr 



2/)2: 



1 - 6ytB cos [(1 - 6) e]g^^ = (1 - 5yeB 



(74) 



The constant LRL vector of the system is obtained simply by substituting its magnitude 
and expression (155|1 for p^ in Eq. ( 15UI) , 



K 



[l-6yiBsm9 

Pr 



2/)2; 



pxe + {i-6ytB 



cos 6 i sin 6 



r'^Pr 



To transform this expression to something of a more familiar form, similar to the classical 
one, it is convenient to introduce, as before, the virtual angle ip = {1 — 6) 9 so that (r, ip) is 
the polar coordinate frame rotating with the orbit. The LRL vector which is fixed in the 
rotating frame is obtained by replacing 6 with ip, 

K' = i^-^f^Bsu.ip^^ ^^ ^^ _ ^^,^,^ ^cos^ _ l^^ ^ 



Pr 



r'^Pr 



which yields, using Eqs. ( 17U]) and (17H]) . 
K' = {l-6)pxi + 



,1 l-a\ E' 
1+ -- 



Ai Mo J c 



^-11 



5a;/i\ K 



^f 



3MaJ 2c2r2 2MocV3 



Using the identity U (ip) K' = U (6) K = Kf which follows from Eq.f l29p . the constant LRL 
vector K is then obtained by rotating K' back to the fixed system. 



K= U (-6) U (^) K' = U{ip-e) K' = U {-59) K' ^ (\ - dOix) K' 



i-5eix] ■ {{i-5)px'i+ 



1 + 



1 l-a\ E' 



H Mo 



r 



5afi\ K^ k(? 

sWo) 2c2r2 ~ 2Moc2r3 



(m 



Finally, the self PB of K are very simply computed by substituting the squared magnitude 



4/)4d2 



K' = {1- 5ytB 



(2a + 3) /ifi;2 
3Moc2 



2,£'+(l-^ 



3fi\ E 



?I2 02 



2 ,,2 



+ yU K 
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in Eq.(l5]) to obtain 

d(KA r / q„\ F'2 



^"^'^I'-AUc^ 



.^jfc^fc ^ (7g) 



thus avoiding much tedious work that would be required for the computation of the PB 
directly from the explicit expression f lTS]) . It is interesting to note that the rhs of Eq. fl76p 
depends only on the energy, and is independent of any other detail of the interaction. As 
well, it is easy to verify that the results of both Subsections coincide for the post-Newtonian 
approximation of a Coulomb system (one charge with infinite mass) . 

VI. THE RELATIVISTIC ORIGIN OF LRL SYMMETRY IN GENERAL 
NEWTONIAN CENTRALLY-SYMMETRIC SYSTEMS 

The last aspect of the LRL symmetry in general rotationally symmetric systems to be 
discussed here is its relativistic origin. After having been shown by Dahl2^ that the classical 
Newtonian LRL vector iQ emerges naturally from the computation of the Lorentz boost in 
post-Newtonian extensions of the Kepler- Coulomb interaction, we show in the following that 
the LRL vector f l43|) of general Newtonian centrally-symmetric systems emerges, in a simi- 
larly natural way, from the Lorentz boost in the corresponding post-Newtonian extensions. 
Consequently, Dahl's original result is not particular to classical Kepler-Coulomb systems 
only, but (at least) to all centrally symmetric ones. 

It should be emphasized that unlike SecJVBj where we looked for the post-Newtonian 
LRL vector, here we look for the Newtonian LRL vector which is derived from the post- 
Newtonian Lorentz boost, in a completely different procedure. 

We start by reviewing the post-Newtonian extensions of Newtonian systems with general 
central interactions^. Consider a 2-particle system with masses mi,m2, spatial coordinates 
Xi,X2; linear momenta pi,P2 and Newtonian central potential Uo{r) with f = xi — X2- In 
post-Newtonian extensions the total linear and angular momenta maintain their Newtonian 
form. 



P = '^Pa=Pl+P2 , J= y^^XgX Pa 



while a scalar interaction Ui {r,Pa) of order 1/c^ is added to the total energy together with 
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the kinetic terms of the same order, 



n ^ 



Pa 



4 
Pa 



+ Uo{r) + -Ui{f,Pa) 



The most general form for the Lorentz boost in the post-Newtonian extension is 



^-E 



"^^+2^ + ^^°^^) 



fa + — ^ - Pt 



(77) 



(78) 



with \1/ {f,Pa) another unknown vector interaction term. 

The properties of the unknowns - f/i and \l/ - are determined by the requirement that P, 
P", J and A^ satisfy, to order 1/c^, the Lorentz-Poincare Lie-Poisson brackets. In particular, 
the PB \n,P°] = P imply 



ES^ + ^^'^n 



dVr, 



—U'^ (r) [f- (vi + V2)] f- -Uo (r) {vi + V2) 



(79) 



with H^ being the Newtonian Hamiltonian (!36|) with potential Uo{r). We need not go 
further into the details of ^ and Ui {f,v), because Eg. (1791) is all that is required for the 
following. 

To obtain the LRL vector, we start from the fact that the Newtonian centre-of-mass 

rriiXi + 1712X2 



Xk 



Mn 



(80) 



{Mo = mi +1712) is constant only in the non-relativistic limit, and look for an internal vector 
R which satisfies, in the centre-of-mass (CM) frame, 

dR _ dX^ _ miVi + m2V2 , . 

Ht ~ dt ~ Wo ^ ^ 

The CM reference frame is defined by P = P1+P2 = 0, without fixing the origin. Substituting 
the post-Newtonian velocities 



gpo 



1 



Pa \ P 



^ + 



1 dUi 



2m?„c^ ) m„ (? 



(82) 



with the the CM condition p\ = —p2 = p and Eq. fj79|) . Eq. flHTj) becomes 
dR {mi — 1712) p'^ ^ ^ 1 -sr^ dUi 

772 2^ 



dt 2m\m\c^ MoC 



dva 



{mi — m,2)p'^ ^ 1 
-p + 



o 2 2 2 ^ ' /i# 2 1 ^^oW[^-(^i + ^2)]r--t/o(r)(t;i + t;2) 

2m,^m,2C'^ MoC^ y2r 2 



MoC^ 



^.H. 



rrii — m,2 
2mim2MoC^ 



V 

- + Uo{r) 
. 1^ 



P--Uoir){r-p)rj-j^^- 



(83) 
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with n = 77111112/ Mo the Newtonian reduced mass. 

Since the rhs is aheady of order 1/c^, we look for solutions of this equation in which 
R is an internal vector of order 1/c^. Newtonian relations and equations of motion are 
then sufficient in the following. As in Dahl's result, Eq. fl83|) is now expected to have two 
independent solutions, with the difference between them being proportional to the LRL 
vector. One solution is obtained from Xn and the Lorentz boost A^. Expressing the particles' 
coordinates in a way similar to the Newtonian relations 









(84) 



the post-Newtonian Lorentz boost fjTHj) becomes, in the CM frame 



^=Y. 



ma + 



p^ 1 



277laC^ 2c^ 



+ IT^Uo (r) 



Xa + ^^(f,p) 



MXn + 



77l2 - 77li p^2 
2fxMoC^ 



[p^ + fiUo (r)] f+ — \I^ {f,p) 



where 



M = Mo + ^ + \Uo (r) 



(85) 



(86) 



is the total relativistic mass. Then the first solution in the post-Newtonian approximation 

is simply 

A^ rrii — 7712 



R\ — Xn — 77 



[p2 + /if/o 



r r 



1 
M„c2 



"^ {f,p) 



(87) 



M 2/iM2c2 

To obtain a second, independent solution for Eq. (l83p . we first employ the Newtonian 
equations of motion and write the equation in the form 



dR 7711—7712 ( d ..^ , _, rrr / N m / m ^ 



dt 2/iM2c2 \dt 



1 dm 

M^~dt 



-n [i^ ■ P) P + 1^ i^Uo)' f] - pr [rUo (r)]" : 



2/iM2c2 [dt 



1 rf^ 

MnC'^ dt 



(88) 



Then, with the vector W {r,p) satisfying Eq. (l42|) . the second solution R2 is identified as 



R2 



mi — 7n2 



2fiMy 

The difference between the two solutions 

7711 — f^2 



r ■ p) p + fi (rUo) f—W 



1 



MoC^ 



r* 



(89) 



i?l — R2 



2/iM2c2 



pxi — firU'^ {r)f+W 



(90) 
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is clearly recognized as being proportional to the LRL vector of the corresponding Newtonian 
system P3|) . We notice that it is independent of the post-Newtonian corrections f/i and ^, 
depending only on the Newtonian limit, so it is the same for all possible post-Newtonian 
extensions of the same Newtonian potential. This result verifies the relativistic origin of the 
LRL symmetry for all Newtonian centrally-symmetric 2-body systems. 

An extra benefit of Dahl's procedure is fixing the preferred direction of the LRL vector. 
Although the classical LRL vector (Q, as historically constructed and used over the years, 
is directed towards the perihelion, there is nothing in the LRL symmetry by itself, as was 
already discussed above (Sec JIIip . that distinguishes any particular direction for it. Choosing 
the direction of the LRL vector in the explicit construction in Sec JIVI towards a closest 
approach was based more on aesthetical reasons (regularity of the coefficients) rather than 
on more profound ones. It is really the relativistic consideration, even at the post-Newtonian 
level, that distinguishes this direction as the preferred one, because this is the direction of 
the vector that appears in Eq. flQOl) . It may be argued, of course, that as an integral, the 
vector i?2 is anyway defined up to an arbitrary addition, so it may result in any desired 
direction. Still, this is an arbitrary addition, while without it the natural direction that 
appears is towards a (generalized) perihelion. 

VII. DISCUSSION AND CONCLUDING REMARKS 

In the present paper we brought and discussed, from several angles, new evidence that 
support the generality and universality of the Laplace-Runge-Lenz symmetry. 
The main aspects that were discussed are : 

1. The independence of the symmetry upon the type of interaction, requiring only internal 
rotational symmetry. This was verified by the propositions of SecUll leading to Eq.(l5]) 
and culminating in Eq.f jTSj) . demonstrating that the symmetry is always o(4) or o(3, 1) 
depending only on the energy state (bound or unbound systems). The dependence of 
the symmetry only upon the energy state is emphasized by the results of SecJV] for 
different types of relativistic systems, where the PB ( l53l) and ( 1761) are independent of 
other details of the orbit (like the status of the potential barrier in Sec. IV AO . except 
the energy. 

26 



These two propositions contain the fundamental ingredients of the internal symmetry, 
incorporating the core property for the generalization of the LRL symmetry to arbi- 
trary rotationally symmetric systems. Eq. i^ thus becomes the fundamental equation 
of the LRL symmetry. 

The generality of the symmetry is strengthened by the fact that the proof of the 
propositions uses only generic properties of the Poisson brackets, independent of any 
particular recipe for their computation. 

2. Definition of the LRL vector - the generator of the LRL symmetry - as a constant 
vector (not only piece-wise conserved) even if the orbit is open, for all centrally- 
symmetric 2-body systems. As an application, the LRL vector in relativistic Coulomb 
systems and post-Newtonian electromagnetic or gravitational systems was computed. 

3. Demonstration of the relativistic origin of the LRL vector in general rotationally- 
symmetric systems. This is a generalization of Dahl's result, which referred only to 
Newtonian 1/r potentials, into all (arbitrary) centrally symmetric potentials. 

Let us discuss now the picture that these aspects combine together. 

The relativistic origin of the LRL symmetry leads us to focus with special interest at its 
appearance in relativistic systems. Relativistic systems, even simpler ones like the relativistic 
Coulomb system or the post-Newtonian systems discussed in SecJVl are characteristically 
endowed with bound states with open orbits. Still, the LRL vector that emerges from 
Dahl's procedure is constant. Thus it was important to insure that we are equipped with 
valid definition for constant LRL vectors, as was done in Sec JIVI 

It has also been shown that although the LRL symmetry, by itself, does not imply any 
particular direction for the LRL vector, the extension to relativistic systems yields a preferred 
direction - towards the perihelia of the orbits, as in the classical Kepler- Coulomb case. This 
follows, for centrally symmetric systems, from the relation between the post-Newtonian 
vector fl90l) and the Newtonian LRL vector fH3|) . 

The LRL symmetry is an internal symmetry - it affects the relative state of the particles 
in a system, but it does not affect the system as a whole (the global, CM motion, remains un- 
affected). In Dahl's procedure (both in its original forn>22 and in the present generalization), 
the LRL vector is derived from the Lorentz boost. This fact points to the possibility that 
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the LRL symmetry is internally associated with global Lorentz transformations, analogously 
to the way that internal rotations are associated with global rotations. In other words, this 
suggests that the LRL and rotational symmetries are attached together internally in the 
same way that global rotations and Lorentz transformations form together the generalized 
rotations in Minkowski space-time (see the diagram below). This also explains why the 
internal symmetry generated by the internal rotations and the LRL vector is o(4) or o(3, 1) 



global Lorentz 

Global space-time symmetry : + 

rotations transformations 



t t 

internal 
Internal symmetry : + LRL 

rotations 

Yet another way to look at the association between the LRL symmetry and the Lorentz 
transformations is by noting that the transformations generated by the LRL vector change, 
for a given value of the total energy E, the internal angular momentum and thus the internal 
configuration of the system, in excellent analogy with the Lorentz transformations changing 
globally the way the system moves as a whole. 

The LRL symmetry is therefore found to be an integral part of the internal relativistic 
symmetry. A detailed discussion of the role that it plays in the internal symmetry of Lorentz- 
Poincare symmetric systems and its implications on the relativistic centre-of-mass is given 
elsewhere^. 

Finally, it should be noted that although the LRL symmetry is known so far to be found 
only in 2-body systems, it follows from Sections [TTl and IIIII that all that is required for its 
existence is rotational symmetry and the existence of PB. Then, following the preceding 
discussion, the LRL symmetry is an integral part of the internal symmetry of the system. 
For 2-body systems, the symmetry generated by i and the LRL vector is the full internal 
symmetry. For larger, A^(> 3)-body systems, the internal symmetry may be larger, but the 
LRL symmetry is expected to be an integral part of it. Therefore, at least in principle, the 
LRL symmetry may well apply also to larger (A^ > 3) systems, as long as they are endowed 
with rotational symmetry. Indeed, two particular cases of the LRL symmetry in many body 
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systems have already been discussed elsewhere^i^. 
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